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1 TITLE 

2 APPROXIMATION METHOD FOR SERIES EXPANSION OF INPUT 

3 FUNCTION AND SYSTEM USING THE SAME 

4 CLAIM OF PRIORITY 

5 [0001] This application makes reference to, incorporates the same herein, and claims all benefits 

6 accruing under 3 5 U.S.C. §119 from my application APPROXIMATION METHOD FOR SERIES 

7 EXPANSION OF INPUT FUNCTION AND SYSTEM USING THE SAME filed with the Korean 

8 Industrial Property Office on 30 April 2003 and there duly assigned Serial No. 2003-27747. 

9 BACKGROUND OF THE INVENTION 

10 Technical Field 

1 1 [0002] The present invention relates to an approximation method for series expansion of an 

12 input function and a system using the same, and, more particularly, to an approximation method 

1 3 and system for series expansion of an input function, which minimizes approximation errors in the 

14 series expansion of the input function of an arithmetic type. 

15 Related Art 

1 6 [0003] Generally, system design starts from upper level algorithm design, in which analysis and 
n verification of the algorithm is performed, and is followed by lower level hardware design. In 
18 order to closely implement arithmetic functions used in the upper level algorithm design in the 
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l lower level hardware design, corresponding function values should be stored in a memory. In this 

^ 2 case, since function values for every input variable should be reserved and stored, the memory 

3 should have a large capacity, resulting in a heavy burden in hardware design. 

4 [0004] Accordingly, in system design, arithmetic functions for input variables are typically 

5 obtained by series expansion, rather than storing every function value for each input variable. 

6 [0005] However, in practice, a series expansion formula cannot be implemented by a 
? combination of an infinite number of terms to reach a nominal function value. Instead, a series 

8 expansion approximation formula, expanded using a finite number of terms, is typically used to 

9 obtain an approximate function value. In this case, an approximation error, or a truncation error, 
l o which represents a difference between the nominal function value and the approximated function 
n value, occurs inevitably. Such an approximation error often causes critical malfunction or 

12 degradation of the system. A series expansion with as many terms as possible can reduce 

1 3 approximation errors, but it will automatically lead to a system of increased complexity required 

14 to implement such a series expansion. 

15 [0006] Let us now assume that an arithmetic function S(x) is expressed in an infinite series 

16 expansion with terms, /o(x) ? /i(x) 5 / 2 ( x ) as follows: 

17 S(x)= f i f n {x)=f 0 {x)+Mx)+f 2 {x)+ (1) 

n = 0 

18 [0007] For example, familiar trigonometric functions often employed in digital system design 
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can be expanded in Taylor series about x=0. In this case, the terms of the trigonometric functions 
have alternating signs (+, -). Some common trigonometric functions expressed using Taylor series 
are as follows: 

cos(x) = Y(-l)" — = 1- — + - (2) 

^0 2n\ 2! 4! K) 



oo x 2n+\ x 3 x 5 

sin(x)= Y(-l)" = x + (3) 

n=0 ( 2 " + 1 >! 3! 5! 



2h+1 3 v 5 7 

arctam>) = YC-l)"- = x _iL + £ — £_ + ... -1<x<1 (4) 

£J 2n + l 3 5 7 



[0008] However, in a fixed-point digital system, series expansion has a finite number of terms, 
thereby causing an approximation error. For example, if a partial sum, 
N 

s N( x ) = ^(~l) n fn ( x ) is used t0 approximate the function, S(x) = (-1)" f n (x) , the 

n=0 n=0 



approximation error is expressed as E^x) = S(x) - S N (x). Since the absolute value Ie n (x) I of the 
approximation error E N (x) increases with the absolute value k I of the input variable x, more terms 
in the series expansion are required to reduce the approximation error. 
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1 [0009] To reduce the approximation error in an alternating series expansion, the Euler 

2 approximation method is frequently used. In the Euler approximation method, the same series 

3 expansion formula as in the Taylor series is applied up to term (N-l), while an Euler 
* 4 transformation formula is applied for the rest of the terms. For example, the function S(x) 

5 presented in equation (1) can be expressed as follows in equation (5): 

6 ZH)V w =/o-/i+/2-"-^-i + I^[AVy V ] J for N = even (5) 

7 wherein, A is called a forward difference operator and has the following characteristics: 
A/jv = /n+1 ~ /n 

8 ^ 2 fN^fN + 2-2f N+l +f N (6) 

A /n s /n+3 ~ 3/n+2 + ?>f N+rfN > etc - 

9 [0010] The Euler approximation method, however, also has a problem in that the rightmost term 

10 should have a large number of terms to achieve an accurate approximation and it works only for 

1 1 N=even. 

12 [0011] Therefore, there is a need to provide an efficient approximation method for series 

13 expansion of functions and a system therefor, which can minimize an approximation error and be 

1 4 simply implemented in practice. 
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1 SUMMARY OF THE INVENTION 

2 [0012] The present invention provides an approximation method for series expansion of 

3 functions, and a system therefor, in which an approximation error of a nominal function value is 
* 4 efficiently minimized. 

5 [0013] According to an aspect of the present invention, there is provided an approximation 

6 method for a series expansion of an input function with a finite number of terms N to minimize 

7 approximation error, comprising the steps of: expanding the input function in Taylor series up to 

8 an (N- 1 )-th term; expanding the input function in Taylor series up to an N-th term; multiplying a 

9 predetermined weight a and the expanded result up to the (N-l)-th term; and combining the 

10 expanded result up to the (N-l)-th term, multiplied by a, and the expanded result up to the N-th 

1 1 term; and dividing the combined result by (a + 1 ). 

12 [0014] It is preferable that a be greater than 0 and less than or equal to 1 . 

1 3 [0015] According to another aspect of the present invention, there is provided an approximation 

14 method for a series expansion of an input function with a finite number of terms N to minimize 

15 approximation error, comprising the steps of: expanding the input function in Taylor series up to 

16 an (N-l)-th term; multiplying an N-th term of the series expansion of the input function and a 

17 predetermined weight value; and combining the expanded result up to the (N-l)-th term and the 

18 multiplied N-th term to be an approximation function / for the series expansion of the input 

19 function. 

20 [0016] It is preferable that the predetermined weight value be Lzlj for 0 < a < 1 . . 

(a + 1) 
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1 [0017] It is preferable that a, obtained for corresponding respective N, be selected to minimize 

2 a maximum of the absolute approximation error. 

3 [0018] It is preferable that a be obtained by: (a) selecting a minimum input in a given input x 
* 4 area; (b) calculating / for the input with the finite number of terms N; (c) obtaining and storing an 

5 error Ie n , x I by subtracting / from a nominal function value of the input x; (d) determining whether 

6 the input x has reached a maximum value in the given input x area, adding a predetermined 

7 increment £ to x if x has not yet reached the maximum value, and repeating steps from (b), (c) and 

8 (d); (e) if x has reached the maximum value, selecting a maximum error value among all the stored 

9 errors of E N , X for all inputs; and (f) searching a to minimize the maximum error value, and storing 

10 a as the weight value for corresponding N. 

n [0019] According to yet another aspect of the present invention, there is provided an 

12 approximation method for a series expansion of an input function with a finite number of terms 

13 N to minimize an approximation error, comprising the steps of: dividing a whole input area into 

1 4 several predetermined sub-intervals : expanding the input function in Taylor series up to an (N- 1 )- 

15 th term in each of the sub-intervals; multiplying an N-th term of the series expansion of the input 

16 function with a predetermined first weight for inputs on left side from a center of each of the sub- 

1 7 intervals, and with a predetermined second weight for inputs on right side from the center of each 

1 8 of the sub-intervals; and combining the expanded result up to the (N- 1 )-th term, and the multiplied 

19 N-th term with the predetermined first and second weights to be an approximation of the input 

20 function in each of the sub-intervals. 

2 1 [0020] It is preferable that the predetermined first and second weights on the left and right side 
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in each of the sub-intervals be selected to minimize a maximum error between the approximation 
of the input function with the finite number of terms N and a nominal value of the input function 
over all the inputs in a corresponding sub-interval. 

[0021] According to still another aspect of the present invention, there is provided a method to 
compensate a carrier frequency offset in an orthogonal frequency division multiplexing (OFDM) 
system, comprising the steps of: estimating the carrier frequency offset s by using a series 
expansion of an arctangent function arctan(x); using the estimated offset to obtain a phase rotation 
value for a first input sample of k=l, wherein sin(27r£) and cos(2;r£) are series-expanded to 

minimize an approximation error; using a phase rotation value for a previous input sample 
including k=l to obtain a phase rotation value for subsequent input sample; and compensating the 
phase rotation values for all the input samples. 

[0022] It is preferable that the estimated carrier frequency offset s be represented by , 



8 = 



In 



arctarr 



Zmy(-i)y*(L-i)) 



> where Re and Im represent a real part and an 



L 



l = \ 



imaginary part of a complex number, respectively, y(i) represents an i-th received sample, L is a 
fast fourier transformation (FFT) size, and £ is an estimated and normalized carrier frequency 

offset of A/T, where T is OFDM symbol duration. 
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[0023] It is preferable that the phase rotation value for the k-th sample be calculated by 



A&T; 
(2n)\ 



2n 



(-1) 



n 



For k = \, cos(Ag)7;)=£ 



sin^)^" 1 ) 



n 



(2*+l) 

(2h + 1)! 



For £ > 2, cos(kAa)T s ) = cos((k - l)A6)T s + A&T S ) 

= cos((k - l)Ao)T s )cos(Aa>T s ) - sin((A: - l)Ad)T s )sm(AST s ) 
sm(kAa>T s ) = sm((k - l)Ad>T s + Ad)T s ) 

= sin((k - l)Ad>T s )cos(Ad)T s ) + cos((£ - l)Ad)T s )sm(Ad)T s ) 

2 [0024] According to still another aspect of the invention, there is provided an approximation 

3 system for carrying out a series expansion of an input function with a finite number of terms N to 

4 minimize an approximation error, comprising an operational processing unit which expands the 

5 input function in Taylor series up to an (N- 1 )-th term, expands the input function in Taylor series 

6 up to an N-th term, multiplies a predetermined weight a and the expanded result up to the (N- 1 )-th 

7 term, combines the expanded result up to the (N-l)-th term, multiplied by a, and the expanded 

8 result up to the N-th term, and divides the combined result by (a + 1). 

9 [0025] It is preferable that a be greater than 0 and less than or equal to 1 . 

1 0 [0026] It is preferable that a obtained for a corresponding respective N be selected to minimize 

1 1 a maximum approximation error. 

1 2 [0027] According to a further aspect of the invention, there is provide an approximation system 

13 for carrying out a series expansion of an input function with a finite number of terms N to 
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minimize an approximation error, comprising an operational processing unit which expands the 
input function in Taylor series up to an (N-l)-th term, multiplies an N-th term of the series 
expansion function and a predetermined weight value, combines the expanded result up to the (N- 
1 )-th term and the multiplied N-th term to be an approximation function / for the series expansion 
function. 

[0028] It is preferable that the predetermined weight value be (z*j for 0 < a < 1 . . 

[0029] It is preferable that a obtained for a corresponding respective N be selected to minimize 
a maximum approximation error. 

[0030] It is preferable that a be obtained by (a) selecting a minimum input in a given input x 
area; (b) calculating / for the input with finite number of terms N; (c) obtaining and storing an 
error E N x by subtracting / from an nominal function value of the input x; (d) determining whether 
the input x has reached a maximum value in the given input x area, adding a predetermined 
increment £ to x if x has not yet reached the maximum value, and repeating steps (b), (c) and (d); 
(e) if x has reached the maximum value, selecting a maximum error value among all the stored 
errors of E N x for all inputs; and (f) searching a to minimize the maximum error value, and storing 
a as the weight value for a corresponding N. 

[0031] According to another aspect of the invention, there is provided an orthogonal frequency 
division multiplexing (OFDM) system for compensating a carrier frequency offset, comprising: 

an estimator for estimating the carrier frequency offset 8 by using a series expansion of a function 
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arctran(x); a first phase rotation calculator for using the estimated offset to obtain a phase rotation 
value for a first input sample of k=l, wherein sin(2;r£) and cos(2;r£) are series-expanded to 

minimize an approximation error; a second phase rotation calculator for using a phase rotation 
value for a previous input sample including k=l to obtain a phase rotation value for subsequent 
input sample; and a compensator for compensating the phase rotation values for all the input 
samples. 

[0032] It is preferable that the estimated carrier frequency offset s be represented by 



. 1 

8 = — arctan 
2n 



llm(y(-i)y*(L-i)) 

i = \ 



lRe(y(-i)y*(L-i)) 

I /=1 



> , where Re and Im represent a real part and an 



imaginary part, respectively, of a complex number, y(i) represents an i-th received sample, L is a 
fast fourier transformation (FFT) size, and s is an estimated and normalized carrier frequency 
offset of A/T. 

[0033] It is preferable that the phase rotation value for the k-th sample be calculated by 
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Ad>T s z 
(2n)\ 



In 



N 



(-1) 



n 



For k = l, cos(A^)=2 



(2n+l) 



sm(AcbT s )=Y J ( ~ l) 



n 



(2n + l)! 



For k>2, cos(kAd)T s ) = cos((* - l)AcoT s + A^) 

= cos((A: - l)Ad)r 5 )cos(A6>7;) - sin((* - l)AaT s )s\n(&®T s ) 
sin(kAd>T s ) = sin((k - l)Ad>T s + A^T;) 

= sin((* - l)Ad)T s )cos(Ad)T s ) + cos((k - l)Ad)T s )sin(Ad)T s ) 



[0034] A more complete appreciation of the invention, and many of the attendant advantages 
thereof, will be readily apparent as the same becomes better understood by reference to the 
following detailed description when considered in conjunction with the accompanying drawings 
in which like reference symbols indicate the same or similar components, wherein: 
[0035] FIG. 1 is a graph showing approximation error characteristics of an arctangent function 
arctan(x), expanded in Taylor series with a finite number N of terms; 

[0036] FIG. 2 is a schematic block diagram for implementing an approximation method for 
series expansion according to the present invention; 

[0037] FIG. 3 A~3C are graphs showing approximation error characteristics of sin(x), cos(x), 
arctan(x), respectively, according to various values of weight a, according to the present invention; 



BRIEF DESCRIPTION OF THE DRAWINGS 
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1 [0038] FIG. 4 is a flowchart for calculating the weight a according to the present invention. 

2 [0039] FIG. 5A is a graph comparing approximation error characteristic of arctan(x) at N=7, 

3 expanded according to Taylor series expansion, Euler approximation method, and an 
* 4 approximation method of the present invention; 

5 [0040] FIG. 5B is a graph comparing approximation error characteristic of arctan(x) at N=8, 

. 6 expanded according to Taylor series expansion, Euler approximation method, and an 

7 approximation method of the present invention; 

8 [0041] FIG. 6 A shows an approximation error characteristic of function arctan(x), the input area 

9 [-1 , 1 ] of which is divided into 4 intervals, each of which is centered at x= -0.75, x= -0.25, x=0.25 

10 and x=0.75, respectively, and separately expanded according to the present invention in each of 
n the 4 intervals; 

1 2 [0042] FIG. 6B shows two comparable absolute approximation error characteristics of function 

13 arctan(x) according to Taylor series expansion and an approximation method of the present 

14 invention, wherein four terms are used, and wherein inputs in both cases are divided into several 

15 intervals in which the arctangent function is separately expanded; 

16 [0043] FIG 7A-7C are constellation diagrams of the Quadrature Amplitue Modulation- 

17 Orthogonal Frequency Division Multiplexing (QAM-OFDM) for 10 OFDM symbols according 

18 to conventional series expansions and the present invention, respectively. 

19 DETAILED DESCRIPTION OF THE INVENTION 

20 [0044] The present invention provides an approximation method for series expansion of 
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functions and a system therefor. 

[0045] FIG. 1 is a graph showing approximation error characteristics, E N (x), of an arctangent 
function arctan(x) implemented by Taylor series expansion with various finite number of terms 
N. From FIG. 1 , it can be noted that the approximation performance is better for inputs near x=0. 
In FIG. 1 , for example, when N=9 and x=-l , it is shown that the approximation error is about 1 .43 
degree. This error may cause a critical system malfunction depending on the system 
characteristics. 

[0046] From the approximation error characteristics in FIG. 1 , the sign of E N (x) for an even and 
odd number of terms alternates at a given x. This is a common phenomenon for the approximation 
error characteristics of all the alternating series functions, such as cos(x), sin(x) functions, 
presented in equations (2) and (3). Table 1 shows the E N (x) sign characteristics of trigonometric 
functions sin(x), cos(x), and arctan(x). 



<Table 1> 





Function 
name 


sin(x) 


cos(x) 


arctan(x) 


N status 


odd 


Even 


odd 


even 


odd 


Even 


x>0 




+ 




+ 




+ : 


x<0 


+ 






+ 


+ 





The trigonometric functions in Table 1 are convergent when approximated in series 
\E (x)\ 

expansion, so thatO < ,' — } \ < 1 and sgn[E N (x)] = sgn{E N ,(x)}. In other words, the signs of 

\E N -M 



E N (x) for even and odd numbers of terms are opposite to each other. By using these 
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characteristics, a new approximation formula for series expansion may be provided as follows: 



[0047] For a given function, the approximation error of the function corresponding to a partial 
sum of (N-l) expansion terms, i.e., S N _i(x), is greater than one having the approximation error of 
the same function corresponding to sum of N expansion terms, i.e., S N (x), and the signs of the 
approximation errors are opposite to each other. To reduce the approximation error of S N .,(x), a 
weight a(0 <a < 1) can be given to S N _,(x), while a weight of 1 is given to S N (x). In equation (7), 
the new approximation formula represents a combination of the partial sum of (N-l) expansion 
terms multiplied by the weight a and the sum of N expansion terms, averaged by a + 1 . Further, 
based on the general formula of alternating series, S N (x) = S^x) + (-1) N / N ( X )> a new 
approximation formula can be obtained from equation (7) as follows: 

S N (x) = S N _ i (x) + -^-f N (x) (8) 

(a + 1) 

[0048] Consequently, in the new approximation formula, where the finite number of expansion 
terms is assumed to be equal to N, the generally known series expansion formula is used for (N- 1 ) 

(-l) N 

terms, while _i I is multiplied by an N-th term, thereby producing an efficient improved 

(a + 1) 
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series expansion formula that reduces the approximation error. Equation (8) is simpler than 
equation (5) which uses an Euler approximation method, and is applicable not only to functions 
having an alternating series characteristic, but also to functions with general series expansion 
having convergence characteristic. 

[0049] FIG. 2 is a schematic block diagram for implementing the approximation method for 
series expansion as described above. The functions of the block diagram can be implemented in 
an operational processing unit (not shown), such as a Digital Signal Processor (DSP). 
[0050] FIG. 3A-3C are graphs showing approximation errors E N (x) of sin(x), cos(x), arctan(x), 
respectively, according to various values of weight a, when the approximation is performed using 
equation (8). Herein, the number of terms N is set, for example to 3, and the optimal value of a 
is obtained by searching for a value which minimizes the maximum of | E N (x) | , the search being 
performed for 0 < a < 1, as represented in equation (9) as follows: 

min{max[|E N (x)|]} (9) 

[0051] In general, the optimal a depends on the number of terms N and the type of function to 
be expanded in series. 

[0052] FIG. 4 is a flowchart for obtaining a weight a for a given N to minimize the maximum 
approximation error in a predetermined interval of input x. 

[0053] A minimum value x min in the predetermined interval is selected in step 402. 
[0054] By using equation (8), an approximation function / = / NjX (N,x) of a given function is 

Page 15 of 35 



PATENT 
P56978 

1 calculated in step 404. 

2 [0055] An approximation error value E N x is obtained and stored by subtracting / from a nominal 
- 3 function value F corresponding series expansion with an infinite number of terms in step 406. 

„ 4 [0056] If it is determined that present x is the maximum value within the predetermined interval 

5 in step 408, the method proceeds to step 412. If not, an accurately chosen increment £ is added to 

6 x in step 410 and the method goes back to step 404. Steps 404 thru 410 are repeated until x is 

7 equal to or greater than the maximum value x max in the predetermined interval. 

8 [0057] If, in step 408, it is determined that x is equal to or greater than the maximum value x max , 

9 in step 412 a maximum absolute approximation error is selected from among all of the 

10 approximation error values of |E Nx | for all x's, stored in step 406. The selected maximum 

1 1 absolute approximation error may be a function of weight a. Therefore, in step 414, a value of a 

12 to minimize the maximum absolute approximation error value can be calculated. 

13 [0058] In step 416, the value of a calculated in step 414 is stored in a memory of the system 

14 which implements this approximation method. 

15 [0059] When this approximation algorithm, based on equation 8, is implemented, a is set in 

16 advance depending on a pre-determined N set through computer simulation. An approximation 

1 7 error of a function expanded in Taylor series using equation (8), based on the block diagram of the 

18 present invention as shown in FIG. 2, by using the optimal weight value a in equation (8), is 

19 effectively reduced, compared with the approximation error obtained using the conventional 

20 Taylor series with a = 0. 

2 1 [0060] FIG. 5 A~5B are graphs for comparing approximation errors of arctan(x) according to the 

Page 16 of 35 



PATENT 
P56978 

1 Taylor series expansion method, the Euler approximation method, and the approximation method 

2 of the present invention. For fairness in comparison, all these methods use the same order of the 

3 highest terms. 

- 4 [0061] With regard to FIG. 5A, both of the approximation methods, the method according to 

5 the present invention and conventional Taylor series expansion, are expanded up to the 7 th term, 

6 i.e., N=7. The Euler approximation method is obtained using the Taylor series expansion up to 

7 3 rd term, and by using the Euler transformation formula for the 4 th to 7 th terms. In FIG. 5A, the 

8 approximation error characteristics, according to the Euler approximation method, appear to be 

9 better than others for inputs whose absolute value | x | is nearer 1 . However, it is shown that the 

10 approximation error characteristics according to the present invention is much better than the 

1 1 conventional Taylor series expansion and the Euler method in overall input areas. 

12 [0062] Referring to FIG. 5B, in the case of N=5 and A 3 , the approximation error characteristics, 

13 according to the Euler approximation method, shows that the error is so large that desirable 

14 performance may not be expected. On the other hand, the approximation error characteristics 

15 according to the present invention are not much different from those in the case of N=7. 

16 Therefore, it is understood that the approximation method according to the present invention is 

1 7 easy to implement and approximates with superior performance regardless of the number of terms 

18 N. 

19 [0063] If it is required to have a uniform approximation error characteristic over a wide input 

20 range, a method in which the input area is divided into several sub-intervals and the new 

2 1 approximation method for series expansion as described above according to the present invention 
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is applied to a function in each of the sub-intervals, can be used (referred as a sub-interval 
divisional approximation method for short). The approximation error becomes much less with the 
sub-interval divisional approximation method. The sub-interval divisional approximation 
algorithm can prevent exponentially increasing approximation errors as input x becomes more 
distant from 0. 

[0064] FIG. 6A shows an approximation error characteristic of arctan(x), which is expanded by 
the sub-interval divisional approximation method, wherein inputs of arctan(x) are divided into 4 
sub-intervals, each of which is centered at x=0.75, x—0.25, x=0.25 and x=0.75, respectively, and 
arctan(x) in each sub-interval is series-expanded according to the present invention. In this sub- 
interval divisional approximation method, each approximation error characteristic does not 
necessarily take the form of alternating series as in the approximation algorithm, the input of 
which is centered at x=0. However, in each sub-interval, a corresponding arctangent function 
arctan(x) approximated by Taylor series expansion still maintains the convergence characteristic. 
Hence, the approximation error can be reduced in each sub-interval by combining the sum of (N- 1 ) 
terms, expanded by Taylor series expansion, and the N th term multiplied by a weight P, obtained 
by using the same principle to obtain a as well as considering the sign characteristic. The sub- 
interval divisional approximation method, to be applied to corresponding sub-intervals, is 
represented as follows: 

S N (x) = S N . 1 (x) + P/ N (x) (10) 



[0065] For each sub-interval, there may be two values of P used. In other words, since the 
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1 approximation error function may have the characteristic of an alternating series in the right or left 

2 directions with respect to a given center in each of the sub-intervals, each sub-interval should be 
* 3 further divided into 2 sub-intervals, for example, [-1, -0.75], [-0.75, -0.5],. ..[0.75, 1], and the 
- 4 optimal value of p is calculated in each of 8 sub-intervals. The optimal value of p in each of the 

5 8 sub-intervals can be calculated in the same manner as that for obtaining a. The sub-interval 

6 divisional approximation method as described above may be implemented in any processor in the 

7 system having an operational processing function. 

8 [0066] FIG. 6B shows two comparable approximation error characteristics of function arctan(x) 

9 implemented by Taylor series expansion and by the method of the present invention, respectively, 

10 input intervals in both cases being sub-divided, and both methods using the same number of 

1 1 expansion terms N=4. In FIG. 6B, p's for the 8 sub-intervals are obtained as 0.888, 1.10, 0.445, 

12 1.56, 1.58, 0.465, 1.1, 0.892, respectively. Since the sub-interval divisional approximation 

13 approach using the Euler approximation method cannot be possible, the approximation error 
H characteristic by the Euler approximation method for comparison are not shown in FIG. 6B. As 

15 shown in FIG. 6B, if an area [-1, 1] (from -1 to 1) is sub-divided into 4 sub-intervals, the 

1 6 approximation error characteristic according to the present invention becomes much less and more 

17 uniform over the input area than the error characteristic according to the conventional Taylor 

18 method. This performance can be more improved if the input area is divided more narrowly. 

19 [0067] Accordingly, an approximation method and apparatus for a series expansion of functions, 

20 according to the present invention, can reduce the approximation error when expanded with a finite 

21 number of terms and has better applicability than Euler approximation method. 
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[0068] Now there will be described a carrier recovery process in an Orthogonal Frequency 
Division Multiplexing (OFDM) system to which an approximation method of the present invention 
is suitably applied. The carrier recovery process in OFDM system often involves Taylor series 
expansion of trigonometric functions, to which the present invention may be applied. 
[0069] For carrier frequency offset estimation in OFDM system, the following equation (11) 
based on the rear part of OFDM symbol and the guard interval is generally used. 



* 1 * 
8 = — arctan<! 

In 



£lm( 7 H)/(I-0) 



Y^{y{-i)y\L-i)) 



(11) 



where, Re and Im represent the real part and imaginary part, respectively, of a complex number, 
y(i) represents an i-th received sample, L is the fast fourier transformation (FFT) size, and 8 is 

the estimated and normalized carrier frequency offset of AfT . When an arctan(x) in equation (1) 

is approximated by Taylor series expansion, the biggest approximation error occurs at x=±l, 

wherein 8 is 8 = 1/8, ±3/8. If it is assumed that the symbol period T = 3.2p,sec as in IEEE 

802.1 la, the biggest estimation error by Taylor series expansion occurs at A/ = ±39.1 KHz, and 
±1 17.2KHz. In addition, if the series expansion of arctan(x) resulted in an approximation error 
of A 0 (° means degree), this means that the normalized frequency offset of A/360 occurs regardless 
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1 of the estimated frequency offset . Since an OFDM system is so susceptible to the frequency offset 

2 that even a small frequency offset can cause a substantial bit-error-rate (BER) degradation due to 

3 inter-carrier interference (ICI), the frequency offset must be kept to a minimum. 

4 [0070] Such an estimated frequency offset 8 = hfT = — is compensated by using a 

5 trigonometric expansion formula like the following equation (12) to estimate the phase rotation 

6 angle of the k-th sample with sample interval T s , that is kAwT s in radians, like the following 

7 equation (1 3). Equation (1 3) is derived from the approximation principle of the present invention. 

8 In following equations (12) and (13), it is noted that initial estimation frequency offset value (for 

9 k=l ) is preferably estimated to prevent cumulative offsets. 

l o sin(a±P) = sin(a)cos(P)±cos(a)sin(P) 

n cos(a±P) = cos(a)cos(p)±sin(a)sin(P) (12) 
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N , nm &o>T s 2n 



For k = l, cos(Afi)r j ) = f; ( 1} (2n)! 



n=0 

(2«+l) 

// „ Aft??; 

^(A^)=2 (_1) (2« + l)! 

For k>2 9 cos(kAd>T s ) = cos((£ - l)Ad>T s + Ad>T s ) 

= cos((* - l)AST s ) cos(A<£7; ) - sm((k - 1) Ad)T s ) sin( A£7; ) 
sin(*A©7; ) = sin((* - 1) AcbT s + Ad>T s ) 

= sin((Jt - \)Ad)T s ) cos(Ad)T s ) + cos((fc - 1) Aa)T s ) $in(Aa)T s ) 

(13) 



[0071] In equation (13), in order to compensate for the phase rotation for the first sample of k=l , 
the series expansion approximation algorithm is used to obtain cos(Am>7^) and sin(Aw>7^) . 

For samples following the first one, k^2, the results of the trigonometric functions obtained for 
previous sample and for the first sample of k=l are used to obtain corresponding values of 
trigonometric functions expanded as in equation (12). In most OFDM system, pilots are inserted 
in a data area despite the data efficiency loss to prevent phase rotation due to the frequency offset 
estimation error. However, in an OFDM system according to the present invention, since a series 
expansion approximation function of the present invention as described above is used to minimize 
the frequency offset estimation error, it is allowed to eliminate such a data efficiency loss factor. 
[0072] FIG 7A-7C are constellation diagrams of the Quadrature Amplitude Modulation- 
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1 Orthogonal Frequency Division Multiplexing (QAM-OFDM), observed for 1 0 OFDM symbols, 

2 according to a conventional series expansion approximation method and the approximation 

3 method of the present invention, respectively. Referring to FIG. 5 again, when the approximation 
. 4 method of the present invention was used at x= -0.94, the worst approximation error occurred. 

5 FIG. 7A and FIG. 7B have been created by performing Taylor series expansion of a function with 

6 a finite number of terms, N=7 and N=8, and resulting approximation errors 0.65° and 0.52°, 

7 respectively. FIG. 7C results from performing the approximation method according to the present 

8 invention. In FIG. 7C, when x= -0.94, the approximation error | E^x) | is 0.047°. From the FIG. 

9 7C, it is shown that QAM-OFDM constellation is placed in optimal position without phase 

10 rotation. 

1 1 [0073] As described above, a new approximation method for Taylor series expansion according 

12 to the present invention enables a simple approximation to be made, is easy to be implemented, 

1 3 and minimizes the approximation error. 

14 [0074] While this invention has been particularly shown and described with reference to 

15 preferred embodiments thereof, it will be understood by those skilled in the art that various 

1 6 changes in form and detail may be made therein without departing from the spirit and scope of the 

1 7 invention as defined by the appended claims. 
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